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Abstract 

, Some physical applications of the Passare-Tsikh solution of a principal 

quintic equation are discussed. As an example, a quintic equation of state 
is solved in detail. This approach provides analytical approximations for 
several problems admitting until now only numerical solutions. 
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*Q_i' 1 Introduction 

^3 , The quintic equation is not of central importance in physics, but in some cases 

it can play a significant role. The first example is provided by Euler's (1767) 
and Lagrange's (1771) researches on the three body problem, in the context of 
celestial mechanics, where, in a particular situation, the distance between two 
\ pairs of bodies satisfies a quintic equation [T] . In modern physics, the equation 

■ of state satisfied by the order parameter describing a 3D liquid mixture is also 

t"**- ! a quintic one [2] . Recently, a quintic equation of state, for pure substances and 

mixtures, has been proposed, as a refined variant of the well-known van der 
\ Waals equation [3 j. The pressure gradient of a fluid in a magnetorheological 

damper satisfies a quintic equation [3] . Other applications occur in the physics 
of molecules, theory of elasticity, etc. (see the references of [5]). An interesting 
example is provided by the phenomenological theory of phase transitins, when 
a physical system is described by a sixth order Ginzburg-Landau expansion of 
the free energy [5J: 

f = -/u + V + (i) 

5_j 2 4b 



O 
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where u is the order parameter and / - an external field. In the most popular 
cases, u is the polarization or magnetization, and / - the electric or magnetic 
field. The equilibrum condition ^ = gives the equation of state: 

f = au + bu 3 + cu 5 (2) 

However, more appropriate for a physical analysis is the inverse relation, 

u = u(f) (3) 

giving the dependance of the order parameter on the external field. In order to 
pass from (2) to (3), we have to solve a quintic equation. 
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Although the roots of a quintic equation are known from the late 1850s 
onwards, their expressions are very complicated, and consequently of limited 
use for analytic calculations. The goal of this paper is to provide a simpler 
approach to this problem. It is essentially based on a recent result obtained by 
Passare and Tsikh [7]. 

The outline of this paper is the following. In the second section, we present 
the most popular methods of solving the quintic equation. The third one is 
devoted to the presentation of the Passare-Tsikh formula for a root of the prin- 
cipal quintic. This formula allows us to express the root of such an equation 
as a series expansion, without making use of a new, and cumbersome, Bring 
(or Tschirnhaus-Bring) transformation. In the next section, using a modern 
presentation of some old results, a Tschirnhaus transformation is constructed in 
detail (in several other papers, it is just postulated, so the reader has a limited 
understanding of its grounds). The fifth section describes how a solution of 
a quintic equation of state - characterizing a system defined by a sixth order 
Landau expansion - can be obtained. In the last section, the relevance of this 
approach, for ferroelectrics and ferromagnets, is analyzed. The possibility of 
obtaining analytic approximations for a quintic equation of state, representing 
a rafined variant of the well knownvan der Waals equation, is also discussed. 

2 An outline of the most popular methods of 
solving the quintic equation 

The first step in solving a "general quintic" 

x 5 + a^x 4 + a^x 3 + a>2X 2 + a\x + a,Q = (4) 

is to reduce it to a simpler form, which does not contain quartic and cubic terms 
- the so-called "principal quintic" (04 = 03 = 0) . This reduction is obtained 
using a Tschirnhaus transformation; in the next section, we shall introduce such 
a transformation in a constructive way, so we shall not give it, now, a formal 
definition. The next step is to solve this simpler "principal quintic" . There are, 
essentially, two approaches to this problem. 

An approach is based on the "Kiepert algorithm" (1878), reformulated in 
modern terms by King and Canfield ([5], [E]). It consists in transforming the 
"principal quintic" into the so-called "Jacobi sextic", whose roots can be ex- 
pressed in terms of some Jacobi functions (theta functions); the roots of the 
"Jacobi sextic" can be then used in order to calculate the roots of the "princi- 
pal quintic" , through an inverse Tschirnhaus transformation. The disadvantages 
of this method are that the transformations involved are quite complicated and 
even the Jacobi theta functions - quite unfamiliar. 

Another approach is to use a Bring transformation (a more subtle - and more 
cumbersome - Tschirnhaus transformation) in order to transform the "principal 
quintic" into a "trinomic quintic" (04 = 03 = a 2 = 0; the commas intend to 
suggest this specific choice of coeficients, giving the simplest trinomic quintic; 
other choices are 03 = = a\ = 0, etc.) . The roots of the trinomic equation: 

x 5 -x~t = (5) 
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can be expressed as generalized hypergeometric functions, 4F3. The "simplest" 
root is given by ([3], eq. (38)): 

„ ,\ 2 3 4 1 3 5 5 5 4N 
*i=*4Fa(-,-,-, -; -, 3 , 3 ;-0= (6) 

£ 9 i 13 i 17 

= t + t 5 + 10- + 15 • 14— + 20 ■ 19 • 18— + ... 
2! 3! 4! 

where we have used Eisenstein's variant [TU] while writting the coefficients of 
the series expansion. 

The other 4 roots are obtained in terms of simple combinations of 4.F3 func- 
tions ([9j, eqs. (39-42)). The disadvantage of this method is that the Bring 
transformation is quite complicated. 

The root (6) was re-obtained recently by Glasser [IT] and Perelomov j!2) . 
using simple methods. These authors gave also, both for equations of lower 
(< 4) or higher (> 4) order, formulae for the roots, expressed as hypergeometric 
functions. 

A remarkable step forward was made by Passare and Tsikh, who were able 
to express a root of a "principal quintic" as a series expansion in its coefficients, 
avoiding, in this way, a cumbersome Tschirnhaus transformation. The main 
point of our paper is to take advantage of the simplicity of the Passare-Tsikh 
formula [7J. 



3 The Passare-Tsikh solution 

Recently, Passare and Tsikh [7J obtained the expression of a root of the principal 
quintic 



as a series expansion: 
x 5 



Bx 5 + Ax 2 + x + 1 = (7) 



j,k>0 ' 



The domain of convergence of this series is given by the condition: 

5 5 |B| 2 - 4 4 |S| + 108|^| 5 - 27\A\ A + 1600|A||B| - 2250|A| 2 |S| < (9) 

It is easy to check that, for A = 0, the root (8) becomes the "simplest" root of 
the trinomic quintic (5); indeed, 

1 2 3 4 1 3 5.5 5 4 
4F3( 5'5'5'5'2'4'4'4^ > {W) 



E 



(5^)! +4fe 



fc!(4fc + l)! 

k>0 y ' 
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Also, for B — 0, (8) becomes 



which is a root of the quadratic equation 

Ax 2 +x + l = (12) 

The convergence condition gives 1 — A\A\ > 0, equivalent to 1 — AA > (the last 
inequality is obviously fulfilled for A < 0), i.e. the condition for x+ to be a real 
root. 

Let us also mention that Passare and Tsikh gave a simple formula for a 
general trinomic equation: 

a + a m x m + a n x n = 0, n > m > (13) 

which can be written, without any loss of generality, as: 

1 + x m + ax n = (14) 

in the form: 

i>o r(i + i±^)fc! 

where e is the m — th root of —1. For m = 1 and n = 2, (15) becomes (11), the 
root of (12). 

Another particular case, m = 3 and n — 5, corresponds, in the classical 
formulation of the Landau theory of phase transitions, to the situation T = T c . 
Indeed, at the critical temperature, the coefficient a in (1) vanishes, a(T = 
T c ) = 0, and the trinomial equation describes the critical curve, which, in the 
case of a magnetic system, is M — M(H) |t=t c - So, the formula (15) allows an 
exact description of such critical behaviours. 



4 The Tschirnhaus transformation 

In a seminal work, published in 1683, Tschirnhaus proposed a method for solv- 
ing a polynomial equation P n {x) = 0, of degree n, which consists, mainly, in 
transforming it into a simpler one, Q n (x) = 0, where the polynomial Q n (x) has 
a simpler form than P n {x), in the sense that the terms of order n — l,n — 2, 
... are removed. However, in practical cases, Tschirnhaus's transformation is in 
general extremely complicated and difficult to use. 

Recently, Adamchik and Jeffrey |13j have given a very useful formulation 
of Tschirnhaus' work (and, also, of Bring's (1786) and Jerrard's (1852) work), 
putting these contributions in a modern, simple and accessible form. In the rest 
of this section, we shall follow closely their paper, in order to reduce a quintic 
equation having the form (2), to a principal quintic. 

Let us note by Xi, i = 1, ...5, the roots of the equation 

x 5 + a 3 x 3 + a\x + a Q = (16) 
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which is similar to (2) and can be obtained from the general quintic (4) putting 
«4 = «2 = 0. Our goal is to reduce the eq. (16) to a "principal quintic": 

y 5 + b 2 y 2 + b lV + b = (17) 

having the roots y k . Let us define: 

5 

S n = S n (x k ) = J2^k (18) 
fe=i 

These sums can be calculated iteratively, using Newton's formula: 

n-l 

S n — -n 5 _„ — S n -jd5-j (19) 

3 = 1 

with a,j = for j < 0. For our particular case, when = a 2 = 0, 

Si = 0, 5 2 = -2a 3 , S 3 = 0, 5 4 = 2aj - 4a u S 5 = -5a (20) 

Also, 

Sa = —S4d 3 — S 2 ai, Si = —5*503 — ^ao, S$ = —S^a^ — 5401 — S^a (21) 

Sg = — 5703 — S5C11 — 5400? SiO — — S%ci3 — 5g&i — S^clq (22) 

Following Tschirnhaus, we shall suppose that the roots yk of (17) are re- 
lated by the roots x k of (16) through the transformation (quadratic Tschirnhaus 
transformation) : 

y k = x\ + ax k + (3 (23) 
Due to the particular form of the principal quintic: 

Si(Vk) =S 2 (y k ) = (24) 

where S n (yk) are defined similarly to (18). Also, 

S 3 (y k ) = -36 2 , S 4 (y k ) = -4&i, S ( y k ) = -56 (25) 

It is easy to find that: 

5 

Si(y k ) = + ax k +P) = S 2 + aSi + 5/3 (26) 

fe=i 

and 

S 2 (y k ) = S i + (a 2 + 2p)S 2 + a[3 (27) 

Using (20), we can express the coefficients a, (3 as simple algebric functions of 
ai and a 3 : 

a=(^a 3 -— ),/?=|a 3 (28) 
5 03 6 

Now, with (25), we can find b 2 , b\, b as functions of a 3 , a\, a using (28) and 
the following three relations: 

S 3 {y k ) = S 6 + 3aS 5 + 3(a 2 + /3)(5 4 + f3S 2 ) + 5/3 3 (29) 
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Si{y k ) = S S + AaS 7 + (3a 2 + 2/3) S 6 + 4a(a 2 + 3/3) S 5 



(30) 



+ (a 4 + 6/3 2 + 12a 2 /3)S 4 + 2/3 2 (3a 2 + 2/3)S , 2 + 5/3 4 

S 5 (yk) = S w + 5aS 9 + 5(a 4 a 2 + 0)S 8 + I0a{a 2 + 2/3)S 7 + 10/3 2 5 6 

+ a{a A + 20a 2 /3 + 30/3 2 )S* 5 + 5/3(a 4 + 6a 2 + 6a 2 /3 + 2/3 2 )S , 4 (31) 

+5/3 3 (2a 2 +/3)S 2 + 5/3 5 

So, we have been able to reduce cq. (16) to a principal quintic, whose coefficients 
b2,bi,bo are determined as functions of a 3 ,ai,a . If we know the roots yt of 
the "principal quintic", we can easely find the roots of (17), using the inverse 
Tschirnhaus transformation defined by (23). 

With these preparations, we can exemplify our approach, solving explicitely 
a quintic equation with physical relevance - the equation of state for a Ginzburg- 
Landau system with sextic anharmonicities. 



5 The solution of the quintic equation of state 

We shall consider the eq. (2), or: 

u 5 + a 3 u 3 + aiu + a = (32) 
which is a particular form of (4), with: 

a 4 = 0, a 3 = - , a 2 = 0, a± = - , a = - — (33) 
c c c 

We shall obtain its roots, following four steps. 

In the first step, we shall transform (32) into a "principal quintic" : 

z 5 + b 2 z 2 + b lZ + b Q = (34) 

using a Tschirnhaus transformation, as described in the previous section. In the 
second step, we shall rescale eq. (34), 

f!> z 2 + &i z + 1 = (35) 

DO O b 

through the transformation: 

z = —w (36) 



to obtain the form (7), with: 



In the third step, we shall use the Passare-Tsikh solution (8), to obtain a (real) 
root of the "principal quintic" obtained from (35) via the transformation (36), 
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having the coefficients (37), as functions of the physical parameters entering in 
the Landau expansion (1). 

Let us describe in some more detail how we can calculate effectively the roots 
of (32). The inverse of the Tschirnhaus transformation (23) is obtained solving 
the equation: 

x 2 +ax+/3-y = (38) 

(to simplify the notation, the indices of x and y have been dropped out) or, 
with x — cz, 

A z + z + l = 0, A = 5— ,c= ,/3 = -p (39) 

or a 

According to [T3] , any of its roots can be used to construct the inverse transfor- 
mation; coosing z+ as this root, it can be written, of course, in terms of radicals, 
or as a series: 

j>0 J w ' 
Expressing iteratively y 5 as given by 

5 = _A 2 _ 1_ _ 1_ 

V ~ B V B V B 

we can write: 

(/?o + y) j = b j0 + ybji + y 2 b j2 + y 3 b j3 + y% A (41) 
where, for j < 4, the form of the coefficients bj n is trivial; for j = 5, 

&50 - Pi - 1 ; 6 B1 = 5/3 4 - 1 ; 6 52 = 10$ - | ; (42) 

&53 = 10$; &54 = 5/?o 

and for n > 5, they can be obtained from the relation: 

6, i+fe = M k b n (43) 

where 6 TO is a column vector with elements 6 m o, ■■■b m 4 and M - a matrix with 
the following non-zero elements: 

1 A 

M„ n = fa, Mib = M 25 = --; A/35 = — 55 M 2i = M32 = M 54 = Af 45 = 1 

(44) 

So, the solution z + takes the form: 

z+ = S + yB 1 + y 2 B 2 + y 3 B 3 + y 4 B 4 (45) 

It would be convenient to express each of y n , 2 < n < 4, in a form similar to 
(8), i.e. y n = X) a jfe'i with given by simple formulas. However, author's 
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attempts of obtaining such expressions failed. Probably, the first step in such 
an endeavour is to identify the Passare-Tsikh series with an already studied 
function, like (for example) the two-variable hypergeometric function. 

In the fourth step, we can obtain the other roots of the "principal quintic", 
constructing the quartic equation whose roots are the four roots of the "principal 
quintic", different of x&, (8). Using Viete's relations, we easily get from (32): 



aci + x 2 + x 3 + x 4 = -x 5 (47) 

xix 2 + xix 3 + X1X4 + x 2 x 3 + X2X4 + X3X4 = x\ (48) 

X1X2XZ + X1X2X4 + X2X3X4 + X1X3X4 = — ao — x\ (49) 

X1X2X3X4 — ao + X5 (50) 
So, xi, X2, X3, X4 are the roots of the quartic equation: 

x 4 — x^x 3 + x\ + (oq + x\)x + ao + x\ (51) 



6 Final comments and conclusions 

It would be interesting to solve the equation of state (2) or (32), in the case 
of a specific ferroelectric. As already mentioned, the coefficient a\ has a linear 
temperature dependence, ao is proportional to the electric field E, and 03 is a 
parameter; so, the coefficients A, B of the Passare-Tsikh solution are polinomials 
in T and E; each term in the series (8) is a polinomial in T and E. In this way, we 
can obtain an analytical approximation for the exact equation of state. However, 
according to a recent compilation of the coefficients of the Landau expansion of 
the free energy [14J . for uniaxial ferroelectrics, like Sro.sBi2.2Ta20g, the sixth 
order term is not yet experimentally measured. A similar situation occurs in 
magnetism 15! . However, in the case of fluids, the quintic equation of state 
allows a more accurate description of the substance near the critical and the 
saturation point [S]; until now, only numerical calculations have been done. Our 
solution provides analytic approximations, useful for a better understanding of 
the physical phenomena. It can be also used in order to extend Sanati and 
Saxena's analysis of the Landau theory [16] . including sextic anharmonicities. 
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